Physical realization of the four color problem in quantum systems 
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A multi-component electron model on a lattice is constructed whose ground state exhibits a 
spontaneous ordering which follows the rule of map-coloring used in the solution of the four color 
problem. The number of components is determined by the Euler characteristics of a certain surface 
into which the lattice is embedded. Combining the concept of chromatic polynomials with the 
Heawood-Ringel- Youngs theorem, we derive an index theorem relating the degeneracy of the ground 
state with a hidden topology of the lattice. The system exhibits coloring transition and hidden- 
topological structure transition. The coloring phase exhibits a topological order. 
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Physics and Mathematics are two disciplines in Exact 
Science which are known to nurture each other. There 
are famous problems in Mathematics which are physically 
realizable. An example is how to determine the config- 
uration of N repulsive point particles confined inside a 
spherical shell. It is then natural to ask whether there 
is a possible realization of other famous mathematical 
challenges. Among them we focus here on the celebrated 
four-color problem, formulated by F. Guthrie in 1852 1]. 
This problem remained unsolved for more than hundred 
years until an affirmative answer was proved by Appel 
and Haken Alongside, a plethora of new mathemat- 
ical concepts were introduced, such as a computer-aided 
proof and Non-deterministic Polynomial-time (NP)-hard 
problems. In the framework of chromatic problems in 
Mathematics, coloring is a passive procedure to the map. 
Namely, the colors are artificially assigned to the regions 
by try-and-error so as to satisfy mathematical rules. A 
more focused question is whether there is an active col- 
oring as a phenomenon in Nature, by which we mean a 
spontaneous emergence of an ordering such as molecular 
crystallization, magnetization, and orbital, whose order- 
ing pattern follows the map-coloring rule. 

In this work we give an affirmative answer to this ques- 
tion in terms of an electron model with several compo- 
nents (colors) defined on an arbitrary lattice. The per- 
tinent ground state exhibits spontaneous ordering which 
follows the map-coloring rule. (Although the Potts model 
is known to exhibit a similar state the mechanism is 
trivial due to the anitiferromagnetic exchange interaction 
between classical spin variables. While ours is a quan- 
tum mechanical one.) Employing the concept of chro- 
matic polynomials and combining it with the Heawood- 
Ringel- Youngs theorem P, 3 , we derive an index the- 
orem, relating the degeneracy of the electronic ground 
state with a hidden topology of the random lattice. Con- 
versely, for a given number of components, we propose 
a constructive method for building a lattice on which 
the exact ground state can be realized. To this end we 
modify the Hajos construction T, 's'l of a random graph. 
The ground state exhibits coloring transition where the 



coloring phase is characterized by a topological order. 
This electronic model, hence, provides a meeting point 
for condensed matter physics, manifold embedding and 
topological graph theory. It is also another route to a re- 
alization of a topological order jHSfe, 
in quantum system. 

The model Hamiltonian (see below) is composed of lo- 
cal non-negative operators. The procedure of construct- 
ing its ground state is very similar to the one employed 
in other models 0, . Electronic realization of such 
Hamiltonian is presented in Refs. [l6l ll^ Hsj l . Here we 
consider a multi-component model. The number of com- 
ponents, e.g. spin or orbital projections, is denoted by 

Nr. 



To expose the cornplete Hamiltonian, the "cell con- 
struction" technique 0, 0| is employed. Consider the 
usual finite random lattice sites A, which can be written 
as A = C U ;B where C f^B = %. (The usual regular and 




FIG. 1: (a) Example of sites and cells. The cells, Ci, C2, 
C3, and Ci are composed of two, three, four, and five sites 
respectively. The open and bold circles denote the capital 
i £ C and border u £ B, respectively. A border can belong 
to several different cells. (See C5, Ce, and C7.) Two cells 
can share several different borders. (See Cs and C9.) Du for 
the site shared by C5, Ce, and C7 is shown by broken line, 
(b) The lattice associated with (a). The lines denote possible 
hopping of electron, (c) The dual graph associated with (a). 
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periodic lattices are allowed.) For each site i € C, we 
define a cell Ci = {i} U Bi where Bi is a subset of B with 
l-Bil > 1 [13. (See Fig.lTta).) The cells C, and Q are 
said to be directly connected if they share at least one 
site (i.e., Ci nCj ^ 0). We assume that every site in B is 
contained in at least one cell and, for any pair of cells Ci 
and Cj, there is a pass of directly connected cells from 
Ci to Cj. A given site u ^ B may be shared by several 
cells and the number of cells to which the site u belongs 
is the coordination number of cell, N^. We assume that 
max{Ny) < Nc- As usual, we denote by Cx^cricx cr) the 
fermion annihilation (creation) operator at the site a; G A 
with the component a, and introduce fermion operators, 

a^,<r = >^c^,. + J2 for leC (1) 

ueBi 

where Du is the set of sites i in C such that u G Ci. 
(See Fig. ^a).) The and /J^ g. create the localized 
single-electron states with the component a in the cell 
Ci and {u} U Du, respectively. Here, we assume that 
the real parameters aI^-* and /i^f'' are nonvanishing and 

they satisfy V"^ + = for any pair of i(e C) 

and u(£ B) contained in the same cell. We consider the 
Hamiltonian 

H = EE^l-^«.- + ^E E 

C7=l u^B iec (T,r—l-a>T 

+ E ^'^i,'^'^u,a, (3) 

iec ueBi (T=i 

where rii^a- = c\ „Ci^a is the number operator, and C/ > 
and V > Q. The first term is the hopping term (see 
Fig. n^b)), the second one is the onsite Coulomb repul- 
sion, and the last one is the density correlation. For 
Nc = 2 with a = (T,i) the system describes fermions 
with spin 1/2. To construct a model with / orbitals one 
set Nc = 2f. With the further assumption such that 
A^"'' are positive for any i e C and u Bi, we can also 
consider another Hamiltonian 

^ = E E ^1-^-.- + E E 

cr=lueB iec <y,T=l-<j>T 

+^EE'^'^-"?^«^<^ (4) 

igC (7=1 

where a[ ,, = Eugb, 

In the following we fix the total electron number Nc = 
J2ia''^i-'^ to \C\ and find the exact ground states of H. 
Consider the states of the form 

*(w)-n"Uio)' (5) 

iec 



where {a} represents a configuration of components 
(cri)igc of the a-operators Let {ctg} be a configu- 
ration satisfying that ct; <jj if the cells Ci and Cj are 
directly connected, and denote by Sq the collection of 
these configurations. Note that a site in B may be occu- 
pied by several (not more than Nc) electrons. Using the 
anticommutation relation {a| g., /3„^o-} = for i G C and 
u G B, we find that the state <I>({ctg}) is a zero mode 
of the Hamiltonians (PJ and (@J and satisfies the relation 
H^{{aQ}) = 0. It is evident that the Hamiltonians are 
non-negative (positive semidefinite) and hence, <I>({(Tg}) 
is an exact ground state. It is possible to prove that any 
ground state can be represented as 

J2 ^({^G})*({aG}) (6) 
{<7G}e5G 

with coefficients 0({(Tg}). 

In order to elaborate upon the ground state and rela- 
tion to the coloring problem, let us fix the terminology. 
A dual graph G(A) is associated with the lattice A in 
such a way that a vertex of G{A) is defined by a cell i of 
A. When the cells Ci and Cj are directly connected, an 
edge in G(A) is defined between the vertices i and j. (See 
Figs. da) and (c).) Note carefully that the dual graph so 
defined is different from the usual dual lattice associated 
with the original lattice A. In standard Graph Theory, 
the lattice and its dual lattice have a one-to-one corre- 
spondence. In the present case, the dual graph is defined 
uniquely for a given lattice, but the inverse is not neces- 
sarily true. Due to this property, we have two kinds of 
redundancy as the lattice is constructed from a given dual 
graph. (See Fig. El) (i) The number of sites composing 
a cell is arbitrary, (ii) The coordination number of cell, 
Nu, is arbitrary even when the coordination number of 
vertices is fixed. This plays a crucial role in the modifi- 
cation of the Hajos construction. As we argue below, the 
existence problem of the ground state is now transformed 
into that of vertex coloring of the dual graph. 

For the case where the dual graph is planar, we can 
regard the cell Ci as the country and the component ai 
as its color. (In this case, the site i is the capital and the 




FIG. 2: Two types of a configuration of cells (a) and (c), 
whose dual graph is the same (b). Then, for a given dual 
graph (b), we have two types of local configuration of the 
cells. The coordination number of the cell, Nu, of the site 
which is shared by Ci, C2, and C3 in (a) is 3. The number of 
the site which is shared by Ci, C2 in (c) is 2. 
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sites in Bi shared by several cells are the border lines.) 
The rule for generating configurations {<tg} is equiva- 
lent to the map-coloring rule. Therefore we call {(Jq} 
as coloring. The answer to the question of whether a 
coloring {ctg} does exist depends on Nc and the struc- 
ture of the lattice. Finding coloring for a given lattice 
is recognized as one of the NP-hard problems. Fortu- 
nately, in the present coloring problem, there are a cou- 
ple of routes to approach it. The first one tells us how 
to look for Nf, for a given lattice. Here the four color 
and Heawood-Ringel- Youngs theorems are applicable re- 
spectively for the two-dimensional planar graph and for 
higher-dimensional random graphs. The second route 
guides us how to construct the lattice for a given Nc- 
Here we need to modify the Hajos construction T, ^ of 
random lattice. 

As a first step, consider a lattice A whose dual graph 
is planar. Employing the four color theorem any 
random planar graph is at least four colorable, and there 
exists at least one ground state of the model with Nc > 4. 
For the model with Nf. < 4, existence of a ground state 
depends on the structure of the lattice. For example, 
the ground state on a lattice whose dual graph is bi- 
partite can be constructed for iVc > 2 because the bi- 
partite graph is two-colorable. When a dual graph is Nc- 
colorable but not (iVc — l)-colorable, we refer to Nc as the 
chromatic number for constructing the ground state and 
denote it by N*. (Note that if a graph is iVc-colorable, 
it is also colorable for the color > Nc Q.) 

Consider now a general random lattice. The problem 
of finding iSg is then equivalent to coloring the random 
hyper-polyhedra. The dual graph G is in general not pla- 
nar and the four color theorem does not apply. Arbitrary 
dual graphs can be embedded into a higher surface with a 
proper genus g Q . They are classified into the torus with 
the Euler characteristic x — 2 ~ 2g, the projective plane 
with X = 1 — 2^, and the Klein bottle with x = — 2(7 [20l |. 
We respectively denote them by T^, , and . Now 
recall the Heawood-Ringel- Youngs theorem 0, 01 . Let 
G be a dual graph on a higher surface S of Euler charac- 
teristic X: and let each cell of the lattice have at most m 
disjoint connected parts. Except for (r7i,x) = (lj2), an 
upper bound on TV* is given by 

_ Tern + 1 + ^(6m + 1)2 - 24x1 



where [ ] is the Gaussian symbol. 

Conversely, for a given Nc, the lattice with N* = Nc 
is obtained by proper recursive use of the following pro- 
cedures (i), (ii) and (iii): (i) Prepare the complete Nc- 
graph, consisting of Nc vertices and edges connecting any 
two of vertices. Replace each vertex by an arbitrary cell 
and identify their sites in such a way that the lattice so 
constructed has the complete A'c-graph as its dual graph, 
(ii) In a dual graph of an already obtained lattice, iden- 
tify two vertices not jointed by an edge and form a new 
dual graph, (ii-a) Replace two cells associated with the 



two vertices identified above by a proper (arbitrary) cell, 
and identify sites so that the resulting lattice maintains 
the connectivity of the new dual graph, (ii-b) If the co- 
ordination number of cell, Nu, of the sites exceeds Nc, 
reduce it by the relation shown in Fig[2Ia)^(b)— >(c). If 
any such reduction attempt fails, increase the number of 
sites in the replaced cell, or go back to (ii) and choose an- 
other pair of vertices for identification, (iii) (modification 
of the Hajos construction): For disjoint dual graphs Gi 
and G2 of lattices which are already obtained, remove an 
edge between the vertices xi and yi in Gi (which means 
the separation of two cells Cx-^ and Cy^^), and remove an 
edge between X2 and y2 in G2. Identify xi and X2, and 
add the edge between yi and j/2, (which means an identi- 
fication of site in the cells Cy-^ and Cy^ ) . Create a lattice 
by following the same ways as (ii-a) and (ii-b). 

Generically, the ground state is degenerate. The degree 
of degeneracy, F, is equivalent to that of the different col- 
oring possibilities of the dual graph G(A). In the coloring 
problem, the number of the different coloring possibilities 
of a graph G is described by the chromatic polynomial 
P{G,x) where x is the number of colors P|. Therefore, 
we obtain the number of zero modes, 

r = P(G(A),A^,). (8) 

For some categories of graphs, the chromatic polynomi- 
als are known. For a general graph, however, derivation 
of its chromatic polynomial is an NP-hard problem. Our 
analysis then shows that the degree of degeneracy is de- 
termined through lO and ((HJ. This index theorem is 
quite intriguing, because the algebraic and topological 
structures of the higher surface determine the properties 
of the interacting electron model. 

As an example, we demonstrate the above procedures 
within a simple lattice. (See Fig. |3Ja).) The dual graph 
is shown in Fig.|3Jb) which is termed as K^,, the complete 
graph of degree five. is irreducibly embedded into 
whose Euler characteristics is x = 2 — 2x1 = 0. From 
eq. we obtain Nc = 7, which is an upper bound of 
N* . The chromatic polynomial for Kr^ is known Q as 

5 

P{K^,x)^\{{x-p+l). (9) 
p=i 

The N* is given by the minimal integer which satisfies 
P{K^,x) > 0, that is, x — 5. Hence, for a given Nc > 5 
the number of zero modes is given by F = Y[p=ii^c''P + 
1). An extension of this procedure to the other lattices 
is straightforward. 

We now discuss two types of exotic phase transitions 
which this system exhibits. Consider an arbitrary lattice 
with N* . Denote by A4 the manifold into which the dual 
graph is embedded irreducibly. 

The first type is a coloring transition '2l'| (For a con- 
text of the Potts model, see ,2^,2^,2^,2^,26||.) Assume 
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FIG. 3: (a) Lattice composed of five cells. Only the borders 
which are shared by more than two cells are shown, (b) Dual 
graph associated with the lattice, (c) Irreducible embedding 
of the dual graph into the two-dimensional manifold. 

that the number of components is controllable by injec- 
tion and absorption. For example, when the components 
are related to electron spin, one can control the num- 
ber of components by applying a magnetic field. And 
when the components are related to an orbital degree 
of freedom, it is controlled by a lattice distorsion. For 
Nc > N*, the ground state is given by ^({ctg}) and the 
system is in a coloring phase. If we reduce the num- 
ber of components and obtain a system with (Nc — 1) 
components on the same lattice, the ground state is also 
described by <i>({crG}) with (Nc — l) components. Within 
this operation, the degeneracy of the ground state under- 
goes a macroscopic change. The corresponding change of 
the system entropy is accompanied by an emission or 
absorption of latent heat. When the number of compo- 
nents is reduced further, the coloring transition occurs 
at Nc — N* because for Nc < N* the ground state is 



not described by <I>({crG})- As the system goes through 
N* the pattern of orbital ordering should undergo a dra- 
matic change. (Mathematically, we define the transition 
point by the smallest zero of the chromatic polynomial 
associated with the dual graph.) The coloring phase, 
Nc > N* , is characterized by a global consistency of col- 
oring. In this sense, the ground state has a topological 
order. 

The second type is a hiddcn-topological structure tran- 
sition. This is a novel structural transition realized for a 
fixed number of components and lattice sites A. In this 
transition there is a variation of the manifold M. into an- 
other manifold M' which undergoes a modification of the 
hopping network. The transition is characterized by the 
set of manifolds [M^M') and, as before, is accompanied 
by a transfer of latent heat. For a given random lattice, 
the ground state properties are essentially determined by 
the topological structure of the higher surface into which 
the dual graph is embedded irreducibly, and not by the 
original lattice structure in real space (hence the adjec- 
tive "hidden" ) . The ground state is then classified by the 
hidden topology. 

Future research directions should include detailed 
studies pertaining to numerous concrete lattices and their 
transitions. The case of infinite lattices is of special inter- 
est in view of the thermodynamic limit and its relation 
to Ramsey theories. 
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